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In supersymmetric extensions of the Standard Model of particle physics (SM), goldstino super- 
partners — scalar and pseudoscalar sgoldstinos — can be light enough for emerging in decays of SM 
particles. Sgoldstino interaction with SM fields is suppressed by the scale of supersymmetry break- 
ing in the whole theory. Hence, searches for sgoldstinos give an opportunity to probe the underlying 
mechanism of supersymmetry breaking. Sgoldstino couplings to SM fields are proportional to the 
supersymmetry breaking parameters — MSSM soft terms — and therefore can lead to flavor violat- 
ing processes in quark and lepton sectors. We consider flavor violating processes involving sgoldstino 
pair production which are driven by sgoldstino couplings proportional to squark and slepton soft 
mass terms, rfiLL ^^d mfj^. We find that present limits on off-diagonal entries in squark and slep- 
ton squared mass matrices allow t-, b-, c-quark and r- lepton decays at levels available for study 
with existing data (BaBar, Belle, CLEOc) and in ongoing experiments (LHCb, CMS, ATLAS). In 
particular, we obtain the following branching ratios Hi {t cSP) < 10~^, Br (r — >■ fiSP) < 10~^, 

Br(B, ^ SP) < 10"*, Br (^B K^*\SP^ < 10~*, Br (D ^ SP) < 10"^ with sgoldstino subse- 
quent decays into kinematically allowed pairs of SM particles 77, e'^e~ , iJ,'^n~, etc. Remarkably, 
the prominent signature of sgoldstino pair production is two muon pairs with pair momenta peaked 
at sgoldstino masses. 

PACS numbers: 12.60.Jv, 13.20.-v, 13.35.-r 



I. INTRODUCTION 

Low energy supersymmetry provides a technically nat- 
ural solution to the gauge hierarchy problem and is still 
very attractive in spite of absence of a clear experimental 
evidence for superpartners. Supersymmetry if it exists, 
must be spontaneously broken and hence there is a spe- 
cial Goldstone supermultiplet, which includes scalar S 
and pseudoscalar P sgoldstinos, goldstino and auxiliary 
bosonic field whose vacuum expectation value F breaks 
supersymmetry in the whole theory. Sgoldstinos being 
massless at tree-level, get masses from higher order cor- 
rections. They can be light if supersymmetry breaking 
happens at relatively low energy, as in models with gauge 
mediation [t] and in no-scale supergravity 

Sgoldstino phenomenology was extensively studied in 
literature (for a brief review see Q) and has got special 
attention after HyperCP anomaly: observation of three 
muon pairs with invariant mass 214.3 ± 0.5 MeV in hy- 
peron decay E — > p/i^/x" The anomaly may be ex- 
plained as pseudoscalar sgoldstino production in E — )■ pP 
with subsequent decay P — 0, [H (scalar sgold- 
stino should be somewhat heavier in this case, ms ^ 
300 MeV, to avoid limits from K — >■ Tr/i^/i"). Sgold- 
stino mass scale lower than the SM superpartner mass 
scale implies that selfinteraction in goldstino sector is 
somewhat weaker as compared to coupling responsible 
for mediation of supersymmetry breaking to the SM sec- 
tor. Sgoldstino explanation can be tested in kaon fS\ and 
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D- and S-meson decays [6| and recent searches with neg- 
ative results 0,13 close some part of the model parameter 
space. 

In this Letter we proceed further with light sgold- 
stino phenomenology concentrating on flavor violating 
processes with sgoldstino pair production. 

Indeed, there are two types of low energy interactions 
between sgoldstino and matter (quark and lepton) fields. 
The type-I involves single sgoldstino field, the type-II uti- 
lizes both scalar and pseudoscalar sgoldstinos. The type- 
I couplings are proportional to so-called left-right soft 
terms in squark (slepton) squared mass matrices, which 
are determined by trilinear soft supersymmetry breaking 
terms. The type-II couplings are proportional to left-left 
and right-right terms of the same matrices. The left- 
right terms are naturally suppressed by fermion (quark 
and lepton) masses with respect to other terms. How- 
ever, the type-II couplings are suppressed by additional 
factor 1/F with respect to type-I couplings. For gen- 
eral low energy supersymmetry breaking models type-II 
couplings might dominate over type-I provided that soft 
trilinear parameters are additionally suppressed. Note 
also, that flavor violating patterns in the two types of 
sgoldstino couplings are generally different. 

Type-I couplings are responsible for sgoldstino decays 
into light SM particles and give rise to single sgoldstino 
production. Possible variants of type-I flavor-violating 
couplings and corresponding phenomenology have been 
studied in detail [1, ■ Type-II flavor- violating terms 
in sgoldstino lagrangian were addressed only in Ref. Q, 
where annihilation of neutral mesons 7r°, Ki^i^s): -C", 
into sgoldstino pair have been studies. In this Letter 
we extend the study of type-II terms to the three-body 
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decays of heavy mesons, t-quark and r-lepton. 



and charged leptons /l., i = 1, 2, 3) read [§| 



II. LAGRANGIAN 

The relevant for present study interaction terms of 
sgoldstino with SM fields (up- and down-type quarks /[/. , 

I 



^ = ^(Sd^P - ((^if +™ff )/l.7^Vl, + « 
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where sum goes over i,j — 1,2,3. Supersymmetry vi- 
olating vacuum expectation value F of auxiliary com- 
ponent of goldstino supermultiplet has the dimension of 
squared mass and the scale of supersymmetry breaking in 
the whole theory is of order ^/F. Parameters m-^j^jj ^y_, 
X = L,R, are MSSM soft supersymmetry breaking terms 
entering squark and slepton squared mass matrices (e.g., 
is the squared soft mass term for superpartner of 



left d-quark, etc). 

There are bounds on off-diagonal elements of matrices 
''^(DUL) ^'^'^ '"^fouL) coming from absence of FCNC 
processes and non-SM contributions to neutral pseu- 
doscalar meson mixings. For numerical estimates of the 
decay rates we use values of soft supersymmetry breaking 
terms which are in agreement with present constraints 
found in literature |13 - [l5l |. It is convenient to define 
the dimensionless variables which parametrize the rela- 
tive strength of flavor violation as 
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further, below we plot decay branching ratios for the spe- 
cial choice of supersymmetry parameters 
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which we call the unitarity limit. |16| We comment on the 
behavior of decay rates with change of model parameters 
in SecEB 

Sgoldstinos are even with respect to i?-parity and 
hence unstable. Apart of ([T]), there are other interac- 
tion terms (see e.g. [3, 9]) responsible for the instability. 
Sgoldstinos decay into SM particles and two-body decay 
mode dominates, see Q for details. Generally, sgoldstino 
decays into 77, e^e~ , fj,'^ii~, tt+tt^ and other hadronic 
modes, if kinematically allowed. Which mode dominates 
depends strongly on the model parameters, but heavy 
leptons in final state are more natural, than light ones 
(for discussion see [EQ)- Thus, for sgoldstino pair pro- 
duction processes one can expect two muon pairs as a 
very promising signature to look for. Other interesting 
signatures include photon pair(s) and light mesons. 



where X = L,R and to refers to a common mass scale of 
superpartners. 

For numerical estimates we have to choose certain val- 
ues of supersymmetric parameters. Without any con- 
firmed evidence for the presence of supersymmetry in 
Nature we are free to choose them at will provided 
no contradictions with experiment. Let us fix general 
scale of SUSY particle masses at to = 1000 GeV, and 
put all i?i?-components of flavor violating terms to zero, 
= . All others are taken as follows 



III. DECAY r 
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= 0.3, 
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Here we consider two body decays of heavy neutral 
mesons into sgoldstino pair, i?" — >■ SP, Bg — >■ SP and 
— > SP. Decay rate of heavy pseudoscalar meson V 
with quark content qiqj (i,j run over quark flavors) into 
sgoldstino pair due to interaction ([IJ is given by the for- 
mula (c.f. Eq.(20) in (1) 
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They are at upper bounds of phenomenologically allowed 
ranges; we have not found in literature relevant limits on 
^LL23 ^^'^ ^LLi3- Hereafter we treat all the values (131)-® 
as reference numbers only. To simplify the presentation 



where M-p and f-p are mass and decay constant of pseu- 
doscalar meson V (see Appendix [A] for the numerical val- 
ues adopted). For decay of antimeson V (heavy meson 
with quark content qjqi) the same formula is valid with 
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substitution i j. For numerical estimates we adopt the 
set of parameter values given in Sec. [Ill The dependence 
of branching ratios of these decays on scalar sgoldstino 
mass ms is depicted in Fig.[Tl where we put pseudoscalar 
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FIG. 1. Branching ratios of decays B° SP and Bs -> SP 
(left) and —>■ SP (right) as functions of scalar sgoldstino 
mass ms for pseudoscalar sgoldstino mass fixed at mp — 
214 MeV and other parameters fixed as explained in Sec.HIl 



Sgoldstino mass to mp = 214MeV, as suggested by the 
sgoldstino explanation of HyperCP anomaly [3, [l] ■ 

Remarkably, the signal grows with sgoldstino mass, as 
the corresponding matrix element exhibits the same be- 
havior, cf . ([Ij , ([7]) . At large mass the branching ratio 
decreases because of limited phase space. 



factors (see Appendix as [13] 

{K°\srb\B°) = [P + p,Y U{ml,) + {P -PiY lA^Az) 
Then the squared matrix element of meson decay reads 
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The decay rate is obtained by integrating |A^p over the 
phase space as explained in Appendix |B] For numerical 
estimates we use the set of parameter values chosen in 
Secini the decay branching ratio for mp = 214 MeV as 
a function of mg is presented on left panel in Figl2] 

Similarly to the two-body decays studied in Sec.lIIII 
the obtained branching ratio grows (but rather slowly) 
with ms in intermediate range of masses, because the 
matrix element contains a term with similar behavior, 
see Eq. (g]). 

Let us proceed with an example of heavy meson three- 
body decay into sgoldstino pair and vector meson, V — >■ 
VSP. We consider decay i?° K*'^SP which goes 
due to the same sgoldstino flavor violating couplings as 
Bs — > SP. We adopt the following notations: P,pi,p2,Pz 
are four-momentums of S", K*^ ^ S and P; e is polariza- 
tion 4- vector of outgoing meson, mf^- = {pi + pj ) are the 
Dalitz variables. We use the following hadronic matrix 
elements [l^ 



IV. DECAYS P VSP AND P VSP 



The same off-diagonal entries of squark squared mass 
matrix which determine the heavy meson two-body de- 
cays into sgoldstino pairs described in Sec. Hill lead to 
three-body decays of the same meson V to light pseu- 
doscalar meson V and sgoldstino pair and to light vector 
meson V and sgoldstino pair. 

As an example of pseudoscalar meson in the final 
state, we consider decay — > K^SP which is driven 
by the same flavor violating squark soft terms which 
trigger Bs — >■ SP. We use the following notations: 
P,Pi,P2,P?, are four-momentums of B^ , K'^ , S and P; 
™fj = {Pi are the Dalitz variables. The relevant 

hadronic matrix element is parametrized by two form 



{K*\srb\B) = (P + {P2^pzf g(™^3) , 

(-P + «+ (923) + (P2 + Ps)^ a- (923) } : 

with form factors described in Appendix|Xl The decay 
amplitude can be written as follows 



M(B ^ K* SP) ^ ~i^x{p2 
+ P2{j>2-PzY{K*\-s^^b\B) 



-p:,Y{K*\si^l''b\B) 



where /3i — 



- 

^023 



^23 



and /32 — 



^23 



We 



present the formulas below for the case of real a and 
/3, which corresponds to P-conservation in squark soft 
terms. For the squared matrix element of three-body 
decay one obtains 
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where M stands for the mass of decaying meson. The 
decay rate is obtained by integrating of |A^p over the 
phase space as explained in Appendix |B] For numerical 
estimates we use the set of parameter values introduced 
in Sccini the decay branching ratio for mp ~ 214 MeV is 
depicted on right panel of Figl2]for kinematically allowed 
values of scalar sgoldstino mass ms- 



where overall numerical factor 3 refers to the number of 
colors, 
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FIG. 2. Branching ratios of B° K°SP (left) and B" 
K*°SP (right) as functions of scalar sgoldstino mass ms- 
Other parameters are as in Fig. [T] 

Several comments are in order. The decay rates for 
charged mesons, e.g. — >■ K^SP, are given by the 
same formulas as presented above. The decay rates for 
(and B^) meson are given by the same formulas with 



replacement m;g. 



'23 



^xx-' 



Likewise it is straight- 



forward to write down formulas for vector meson decays 
into sgoldstino and for various three-body decays of pseu- 
doscalar and vector £)-mesons and Be mesons. 



^23 



U23 



for t ^ cSP 



4F2 



for t uSP. 



and 71 = , 72 

and 71 = — "^pi , 72 

For numerical estimates we use the set of values chosen 
in Secini For the mass and width of top-quark we take 
172.0 GeV and 1.3 GeV respectively. Integrating |A^p 
over the phase space as explained in Appendix [B] one 
obtains the decay branching ratios which are presented 
in Fig. [3] (left) as functions of ms for mp = 214 MeV. 
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FIG. 3. Branching ratios for t c{u)SP (left) and r 
fi{e)SP (right) as functions of ms- Other parameters are as 
in Fig. [U 



V. t c{u)SP AND r ^ /i(e)5'P 

Flavor-violating couplings of sgoldstino to top-quark 
in Eq. ([T]) result in decays t — cSP and t — > uSP- 
To describe it below we use the following notations: 
P,pi^P2,P3 are four-momentums of t-quark, lighter up- 
type quark (c or u), S and P; M denotes mass of top- 
quark, mf 



squared matrix element of the three-body decay can be 
written as 



Decays r — >■ ^SP and r eSP can be treated 
in the similar way. Now the notations are as follows: 
P,pi,P2,p-i are four-momentums of r-lepton, light lep- 
ton I (muon /i or electron e), S and P; M stands for 
T-lepton mass, mf^ = {pi+Pj) are the Dalitz vari- 
ables. The squared matrix element of the three-body 
decay is presented by Eq. ^ without color factor 3 and 



{pi+Pj) are the Dalitz variables. The with 71 = 



^^23 ^23 
4_F2 



72 



ip2 



for decay 



T ^SP and 7i = 4^ 2 "'^ , 72 = ^^"^4^2 
for decay r — >■ eSP- The decay rates are obtained by 



l-^l = 3 ■ [(I71P + I72P) Tl + (I72P — I71P) T2] , (9) integrating \M\'^ over the phase space as explained in 
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Appendix |B] For numerical estimates we adopt the pa- 
rameter values given in Sec. [Ill the decay branching ra- 
tios for mp = 214 MeV are presented in Fig. [3] (right) as 
functions of mg. 



VI. SUMMARY 

To conclude we have studied sgoldstino pair produc- 
tion by heavy quark and r-lepton decays. The obtained 
estimates of corresponding branching ratios for B- and 
D-mesons, i-quark and r-lepton decays show that the 
processes are available for study already with collected 
statistics in Belle, BaBar, LHCb, CLEOc. The numeri- 
cal results are presented for a particular set of supersym- 
metry parameters, see Sec. |lTl These are mixing angles 
in squark (slepton) squared mass matrices 6, superpart- 
ner mass scale rh and supersymmetry breaking parame- 
ter F. All the branching ratios discussed above scale as 
c>c ■ rh'^/F'^, which allows to give predictions for other 
phenomenologically viable values of the model parame- 
ters. 

Note that generally, soft supersymmetry breaking 
terms violates CP-symmetry, that, in particular, splits 
rates of sgoldstino pair production in particle (meson) 
and antiparticle (antimeson) decays, e.g. 

Br {Bs SP) ^ Br {B, SP) . 

Similar statement holds for the three-body decays. In 
case of two-body decays the splitting (so-called asymme- 
tries) can be estimated with the help of decay rate for- 
mulas presented in Sec. lIIII Three-body decay rates with 
CP-violation can be calculated with help of formulas for 
amplitudes given in Sec. lIVI 
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Appendix A: Meson Form factors used in numerical 
estimates 

For meson leptonic decay constants we use fo — 
207 MeV fB = 190 MeV [11, /s, - 231 MeV [H. 
For other meson form factors we adopt the universal form 
of dependence on momentum transfer [T7| 



with values of dimensionless parameters listed in the Ta- 



F(G) 



l-aiqyM^) + b{qyM^Y 



F 




pBK 


yBK 


^0 






FiO) 


0.35 


0.35 


0.31 


0.31 


0.26 


0.24 


a 


1.58 


0.71 


1.79 


1.68 


0.93 


1.63 


b 


0.68 


0.04 


1.18 


1.08 


0.19 


0.98 



TABLE I. Numerical coefficients in meson formfactors 117 



mulas in the main text, as follows Uj 



- ml 



V{ql,)^-{M + m,)g{ql,)^ Mil.) ^ -J^^^ 



2m 1 



923 j 



where 



^3 (9I3) = ^1 (923) - ^2 ('?2'3) 
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Appendix B: Three-body piiase space integration 

To calculate the rates T of the three-body decays de- 
scribed in the main text one has to integrate presented 
there squared matrix elements over the three-body 

phase space of corresponding final state. For the decay 
of particle of mass M into three particles of masses rrii 
with 4- momenta pi, i ^ 1,2,3, one has the general for- 
mula (see e.g. Sec. 39 in Ref.[20|) in terms of the Dalitz 
variables to? ■ = {pi + Pj) , 



1 1 



(mi2)7nax 



(™23)" 



dml,\M\' 



(2^)3 32M3 7(,„.^)_ A™i3). 

where {ml2)niin = (™1 + ™2)^ (™?2)ma2; = {M - ma)^. 



and E* = "M2-"M+'"2 ^* ^ 

^ 2mi2 ' 3 
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